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Abstract 

We introduce the total s-energy of a multiagent system with time-dependent hnks. This 
1-H provides a new analytical lens on bidirectional agreement dynamics, which we use to bound 

the convergence rates of dynamical systems for synchronization, flocking, opinion dynamics, and 
social epistemology. 
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1 Introduction 

(N 

We introduce an analytical device for the study of multiagent agreement systems. Consider an 
infinite sequence of graphs, Go,Gi,G2, ■ ■ ■, each one defined on n nodes labeled 1, . . . , n. We assume 

that each graph Gt is embedded in Euclidean d-spacc and we let Xi{t) G R'^ denote the position of 
node i at time t. The total s-energy E{s) of the embedded graph sequence encodes all of the edge 
^ lengths: 

§ E{s) = Y^ Yl \\xi{t)-xj{tm, (1) 

7^. t>o {i,j)eGt 

. where the exponent s is a real (or complex) variable. The definition generalizes both the Dirichlet 

form derived from the graph Laplacian and the Riesz s-energy of points on a sphere. Sometimes, 
^ variants of the total s-energy are more convenient; for example, we will use the kinetic s-energy, 

n 

K{s)=YY.\\^^{t + l)-X,{tm. (2) 
t>0 i=l 

Note that these definitions make no assumptions about the temporal network, which is the name 

given to a graph sequence sharing the same node set. There is no reason to think that s-energies 
should even be finite, let alone useful: for example, -E(O) is usually infinite. In fact, it is immediate 
to embed a temporal network so as to make its total s-energy diverge everywhere, so one clearly 
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needs assumptions on the embeddings. In this paper we consider the case of multiagent agreement 
systems [33], which we define in the next section. There are two kinds: bidirectional and nonbidirec- 
tional. We consider only the former type in this work. We thus assume that each Gt is undirected, 
meaning that if is an edge then so is (The directed case is quite different and warrants 

a separate treatment.) 

We use the total s-energy to bound the convergence rates of classical systems for opinion dynam- 
ics ( ^2.1[ ), social epistemology ( ^2.2[ ), Kuramoto synchronization ( ^2.3[ ), and bird flocking ( ^2.4[ ). We 
deal only with discrete-time dynamical systems or, as in [33j, time-1 maps of continuous systems. 
We also improve a classic bound for the products of stochastic matrices ( ^2.5[ ). 

Our proofs are algorithmic and make no use of algebraic graph theory (with a single exception 
for reversible systems where we turn to a standard £2 argument ^ 1.1 ). In particular, the proofs focus 
on the agents rather than on the matrices governing their dynamics. In fact, the proofs themselves 
can be viewed as dynamical systems which embed convergence measures directly within the inter- 
agent communication. We hope this perhaps opaque comment finds clarification below and that the 
benefits of an algorithmic approach to multiagent dynamics becomes apparent [6l [8] . 



1.1 Multiagent Dynamics 

Moreau ^33j introduced a geometric framework for multiagent agreement dynamics of appealing 
generality. He established convergence criteria based on connectivity in both the directed case 
and the undirected one. We seek to analyze the dynamics of bidirectional systems without any 
connectivity assumptions: specifically, our goal is to provide bounds on convergence rates that hold 
in all cases. 



Bidirectional agreement systems. The one-dimensional case features all of the ideas and it 
is straightforward to extend our analysis to d > 1; we briefly mention how to do that below. For 
simplicity, therefore, we assume that d = 1. The model involves n agents located at the points 
xi{t), . . . , Xn{t) in M at any time t > 0. The input consists of their positions at time t = 0, together 
with an infinite sequence {Gt)t>o of undirected graphs over n > 2 nodes (the agents); each node 
has a self-loop. These graphs represent the various configurations of a communication network 
changing over time. The sequence need not be known ahead of time: in practice, the system will 
often be embedded in a closed loop and the next Gt will be a function of the configuration at time 
0, . . . ,t — 1. The strength of the model is that it makes no assumption about the generation of the 
temporal networks nor about their connectivity properties. In the case of directed graphs, such a 
level of generality precludes blanket statements about convergence; bidirectionality, on the other 
hand, allows such statements. The neighbors of i form the set Ni{t) = {j \ {i,j) £ Gt}, which 
includes i. At time t, each agent i moves anywhere within the interval formed by its neighbors, 
though not too close to the boundary: formally, if mi^t is the minimum of {xj{t) \j G Ni(t)} and 
Mi^t is the maximum, then 

(1 - p)m,^t + pMi^t < Xi{t + 1) < pnii^t + (1 - p)Mi^t, (3) 

where 0<p<l/2is the (time-independent) agreement parameter, fixed once and for all. All 
the agents are updated in parallel at each step t = 0,1,2, etc. We conclude the presentation of 
bidirectional agreement systems with a few remarks. 

• The model describes a nondeterministic dynamical system. This refers to the fact that the 
sequence of graphs, as well as the particular motion of the agents, are left completely arbitrary 
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within the constraints imposed by ([3|: they could be decided ahead of time or, as is more com- 
mon, endogenously in a closed-loop system; we give several examples below. The embedding 
at time is provided as input and, from then on, all subsequent embeddings are generated by 
the system itself in abidance by rule ([s]) . It may sound surprising at first that one can prove 
convergence in the presence of such high nondeterminism and without the slighest assumption 
about connectivity. 

• Bidirectionality does not imply symmetry among neighbors. In fact, the behavior of neigh- 
boring agents may be completely different. The condition p > is essential. Without it, a 
two-agent system with a single edge could see the agents swap places forever without ever 
converging. This simple example shows that one may legally move the two agents toward each 
other so that their distance decreases by a factor of merely 1 — 2/9 at each step. This shows 
that no worst-case convergence rate can be faster than e~'^^^. 

• There are several ways to extend the model to higher dimension. Perhaps the easiest one is 
to assume that agent i is positioned at Xi{t) G M'^ and then enforce ([s]) along each dimension. 
This is equivalent to having d one-dimensional systems sharing the same temporal network; 
it is the method we use in this paper. A different, coordinate-free approach stipulates that 
agent i may move anywhere within the convex hull Cj (t) of its neighbors { Xj [t) \j € Ni (t) } , 
but not too close to the boundary (Figure [T]) . This requires shrinking Ci (t) by a factor of 1 — p 
centrally toward a well-chosen center: for example, the Lowner-John center oi Ci{t), which is 
uniquely defined as the center of the minimum- volume ellipsoid that encloses Ci{t) |14j . 




Figure 1: The agent can move anywhere inside the pentagon but may not touch the thick boundary. 

Much of the previous work on agreement systems has been concerned with conditions for con- 
sensus (ie, for all agents to come together), beginning with the pioneering work of [41^ W2\ and 
then [H El m m [IHl [m [211 [331 IM] . Bounds on the convergence rate have been obtained under vari- 
ous connectivity assumptions |5l [51] and for specialized closed-loop systems [71 [30] . The convergence 
of bidirectional agreement systems can be derived from the techniques in [18^ [271 [33] . Bounding the 
convergence rate, however, has been left open. This is the main focus of this paper. Before stating 
our results in the next section, we discuss a few extensions of the model. 

The fixed-agent agreement model. We can fix one agent if we so desire. By this, we mean 
skipping the update rule at an arbitrary agent iq, selected ahead of time — or equivalently, directing 
all edges incident to iq toward that node. To see why, create the point reflection of the n — 1 mobile 
agents about iq to create a bidirectional system of 2n — 1 agents. Figure [2] illustrates this process 
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in two dimensions for visual clarity. We duplicate each graph Gt, with the exception of the fixed 
agent iq. In this way, at time t, each edge is given a duplicate Placing the origin of a 

Cartesian coordinate system at XiQ{0), we position agent i' at time so that Xj'(O) = — Xi(0), which 
inductively implies that Xj'(t) = —Xi{t) for all t > 0. No edges connect the two copies of the original 
graphs. Every mobile agent (and its reflected copy) mimics the behavior of its counterpart in the 
original n-agent system while respecting ([s]). The fixed agent always lies at the midpoint of the 
smallest interval enclosing its neighbors; therefore, it does not need to move, even for the maximum 
value of p allowed, which is 1/2. To summarize, any n-agent agreement system with one fixed agent 
can be simulated with a (2n — l)-agent bidirectional agreement system with the same value of p and 
at most twice the diameter. We apply this result to truth-seeking systems in §2.2[ 




Figure 2: Reflecting the system about the agent lo that we wish to fix. 



Reversible agreement systems. Assign to each agent i a time-independent motion parameter 
Qi ^ l^j(^)l ^-iid define the mass center of the agent's neighbors as 



A reversible agreement system satisfies the transition: 

Xi{t + 1) = Xi{t) + ^-^^ 



Pi{t) - xi{t)y 



The agents obey the dynamics x{t + 1) = P{t)x{t), where 



Ptjit) 



l-imt)\-l)/qi ifi=j; 

1/qi ii i ^ J G Ni{t); 

else. 



(4) 



A quick examination of Q shows that the dynamics, indeed, defines an agreement system with 
parameter p = 1/ maxj qi. Why reversible? We take note of the identity qiPij{t) = qjPji{t). This is 
the standard balanced condition of a reversible Markov chain, with (qi) in the role of the stationary 
distribution (up to scaling). Indeed, we easily verify that the sum qiXi{t) is independent of t and 
that a lazy random walk in a graph is a special case of a reversible agreement system. The latter is 
much more general, of course, since the graph can change over time. We note that, if each node has 
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its own degree fixed over time, tlien moving each agent to the mass center of its neighbors satisfies 
reversibihty and hence, as we shall see, fast convergence. This is equivalent to setting qi = \Ni{t)\. 

The definition of a reversible system is simple yet somewhat contrived. Does the concept lend 
itself to intuition? Yes. At each step, each agent picks a set of neighbors (the graph Gt) and slides 
toward their mass center — but not all the way! The agent might have to stop before hitting the 
mass center. When? This is where the qi's come in. They ensure that the gj-weighted mass center of 
the whole system stays put. Not to have that mass center wiggle around and produce exponentially 
small coefficients is a key reason why reversibility implies faster convergence. We flesh out this 



intuition in ^ 1.3 



The matrix approach. It is customary to model agreement systems by using products of stochas- 
tic matrices: x{t+l) = P{t)x{t), where x{t) = {xi{t), . . . , Xn{t))'^ and P{t) is a row-stochastic matrix 
whose entries Pij{t) are positive for all i,j with j £ Ni{t). Bidirectionality means that Pij{t) and 
Pji{t) should be both positive or both zero, which is a form of mutual confidence [27j. Typically, 
one also requires a uniform lower bound on each nonzero entry of the matrices. We observe that 
condition ([s]) is not nearly as demanding: all we require is that, if the agent i has at least one 
neighbor (besides itself), then the entries corresponding to the leftmost and rightmost neighbors l{i) 
and r{i) should be at most 1 — p. These conditions have a natural interpretation that we summarize 
below: for all t, 

{Mutual confidence: No pair Pij{t),Pji{t) has exactly one zero; 
No extreme influence: For any nonisolated agent i, max{pjj(j)(t),pjf.(j)(t)} < 1 — p. 

Conditions ([S]) are weaker than the usual set of three constraints associated with the bidirectional 
case [IHllST], which, besides mutual confidence, includes: self-confidence (nonzero diagonal entries) 
and nonvanishing confidence (lower bound on all nonzero entries). Our model requires bounds 
on only two entries per matrix row. Previous work [21 EJ [27] highlighted the importance of self- 
confidence {pii{t) > 0) for the convergence of agreement systems. Our results refine this picture: 
To reach harmony in a group with bidirectional communication, individuals may be influenced ex- 
tremely by non-extreme positions but must be influenced non-extremely by extreme positions (rrii^t 
or Mi^t)- In the case of a two-agent system, this maxim coincides with the need for self-confidence; 
in general, the latter is not needed. We conclude this comment about the matrix representation 
of agreement systems by emphasizing that the total s-energy seeks to move the focus away from 
the matrices themselves and, instead, reason about the agents' motion in phase space and their 
temporal communication network. 



Random walks and ergodicity. At the risk of oversimplifying, one might say that to understand 
agreement systems is to undestand backward products of stochastic matrices, 

P(t)P(t-l)---P(l)P(0), 

as t grows to infinity. Forward products P(0)P(1) • • • P(t), for t — t- c«, are different but much can be 
inferred about them from the backward variety. A forward product of stochastic matrices models 
a random walk in a temporal network: imagine walking randomly in a graph that may change at 
every step. These have been studied by computational complexity theorists, who call them colored 
random walks |1UI [TT] . This connection suggests that a complete theory of agreement systems would 
need to include, as a special case, a theory of discrete-time Markov chains. As we shall see, the total 
s-energy allows us to retrieve classical mixing bounds for random walks in undirected graphs. 
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Figure 3: The total s-energy is a global instrument to track how fast the polytopes shrink. 



A general principle behind the convergence of products of stochastic matrices is that, if all goes 
well, as t grows, the product will tend to a matrix of rank one or a (possibly permuted) block- 
diagonal matrix with blocks of rank one. Many analytical devices have been designed to keep track 
of this evolution, most of which fall into the category of ergodicity coefficients [38]. There is a simple 
geometric interpretation of this which is worth a short detour. From a stochastic matrix such as 
P(0), construct a convex polytope by taking the convex hull, denoted convP(O), of the points formed 
by the rows of -P(O): here, each row forms the n coordinates of a point. When we multiply P{1) 
by -P(O), each row of the product is a convex combination of the rows of P(0), so the corresponding 
point lies inside the convex hull convP(O); therefore conv{P(l)P(0)} C convP(O). A coefficient 
of ergodicity is typically a measure of how quickly the nesting shrinks: it might keep track of the 
width, diameter, volume, or any other "shrinking" indicator: 

conv{P(t) • • • P(0)} C conv{P(t -!)••• P(0)} C • • • C conv{P(l)P(0)} C convP(O). 

For example, when all the matrices are identical (as in a Markov chain), the spectral gap gives us an 
^2-norm tracker of the contraction. What all coefficients of ergodicity have in common is that they 
are local instruments: they measure the contraction from one step to the next. The total s-energy, 
instead, is a global instrument. It monitors the shrinking over all time steps in a global fashion: 
the parameter s plays the role of frequency in Fourier analysis and allows us to choose the correct 
"frequency" at which we want to monitor the shrinking process. This gives us Chernoff-like bounds 
on the distribution of the edge lengths. 

1.2 The Total s-Energy 

There is no obvious reason why the total s-energy, as defined in ([T]), should ever converge, so we 
treat it as a formal series for the time being. We prove that it converges for any real s > and 
we bound its maximum value, En{s), over all moves and n-node graph sequences. We may assume 
that all the agents start out in the unit interval (which, of course, implies that they remain there 
at all times). The justification is that the total s-energy obeys a power-law under scaling: x i— )• Cx 
implies that En{s) i— )• C'^En{s). We also assume throughout the remainder of this paper that p is 
smaller than a suitable constant. All the proofs of the results below are deferred to ^ 
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Theorem 1.1. The maximal total s-energy of an n-agent bidirectional agreement system with unit 
initial diameter satisfies 



En{s) < 



p-0{n) f^^ g^i. 

^i_„^_„2_0(i) for 0<s<l. 



There is a lower bound of 0{p) l-^/^J on En{l) and of s En{s) for n large enough, any 

s < sq, and any fixed sq < 1. 

The asymptotic notation hides the presence of absolute constant factors. For example, p~^^'^\ 
0(p)~l-"/2J^ and p~^^"'^ mean, respectively, at most p""^", at least (6/?)~l-"/2J ^ and at least p'^", for 
some suitable constants a,b,c > 0. Since no edge length exceeds 1, En{s) < En{l) for s > 1, and so 
the theorem proves the convergence of the total s-energy for all s > 0. 

When the temporal network always stays connected, it is useful to redefine the total s-energy 
as the sum of the s-th powers of the diameters. Its maximum value, for unit initial diameter, is 
denoted by 

Enis) = J2 (diam{2;i(t), . . . ,x„(t)})'. 

In dimension d = 1, the diameter is the length of the smallest enclosing interval. The following 
result is the sole breach of our pledge to avoid any connectivity assumption. 



Theorem 1.2. The maximal diameter-based total s-energy of a connected n-agent reversible agree- 
ment system with unit initial diameter satisfies 



n-^E^{s)<E^{s)<—( — ) 

s \ p / 



for all < s < 1. 




Figure 4: By analytic continuation, the maximum total s-energy of a two-agent system is a meromorphic 
function over the whole complex plane; the function is depicted in absolute value, with the real axis stretching 
to the right. 
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We proceed with general remarks about the function E{s). All of the terms in the series are 
nonnegative, so we can assume them rearranged in nonincreasing order. This allows us to express 
the total s-energy as a general Dirichlet series: 



(6) 



k>l 



where Afc = — Indk and is the number of edges of length dk- Thus, E{s) is the Laplace transform 
of a sum of scaled Dirac delta functions centered at x = A^. This implies that the total s-energy 
can be inverted and, hence, provides a lossless encoding of the edge lengths. We show that E{s) 
converges for any real s > 0. By the theory of Dirichlet series [15], it follows that E{s) is uniformly 
convergent over any finite region D of the complex plane within 5ft(s) > r, for any r > 0; furthermore, 
the series defines an analytic function over T>. It is immediate to determine the maximum s-energy of 
a 2-agent system with unit initial diameter. For p = 1/2 — l/2e, E2{s) = ^^^(l — 2/j)** = 1/(1 — e~*); 
therefore, writing s = x + iy, it satisfies (Figure |4]): 



The singularities are the simple poles s = Inik, for all k. The maximal total s-energy can be contin- 
ued meromorphically over the whole complex plane. Note that this is obviously false for nonmaximal 
s-energies: for example, the function g-sfc! jg valid total s-energy, but its singularities form 
a dense subset of its line of convergence (the imaginary axis), hence an impassable barrier for any 
analytic continuation into K(s) < 0. 



Figure 5: ^-Convergence is reached when the agents fall within groups with disjoint enclosing intervals of 
length at most e and no further interaction ever takes place between distinct groups. 

1.3 Convergence 

Bounding the convergence rate of agreement systems faces the obvious difficulty that an adversary 
can always make the temporal network free of any edges joining distinct nodes for as long as it 
wants, and then, at some point far into the future, add all (2) edges permanently to the temporal 
network in order to make all the agents cluster around the same point. How then can one hope to 
bound the convergence time since it can be arbitrarily large yet finite? The total s-energy is meant 
to deal with precisely that problem. Given < e < 1/2, we say that a step t is trivial (where e is 
understood) if all the edges of Gt have length at most e. Recall that Gt always has n self-loopsP] 




The communication count G^ is defined as the total number of nontrivial steps. Intuitively, it is a 
way of ignoring microscopic motions. The system is said to e-converge if the n agents can eventually 
be partitioned into subsets with disjoint enclosing intervals of length at most e and, from that point 
on, no further interaction ever takes place between distinct subsets (Figure [5]). Consensus refers to 
the case of a one-set partition. Visually, e-convergence means that the system eventually freezes. 
We have this obvious relation between e-convergence and triviality: 

self-loop always has zero length but an edge may be of zero length without being a self-loop. 



E2{s)\ = l/Vl-2e-^'cosy-Fe 







1 
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Fact 1.3. An n-agent bidirectional agreement system (n — l)e-converges by the time its last non- 
trivial step has elapsed. If the temporal network remains connected at all times, then the system 
(n — l)e-converges to consensus within Ce time. 



Theorem 1.4. The maximum communication count Ce{n) of any n-agent bidirectional agreement 
system with unit initial diameter satisfies 

0(p)-L"/2j log 1 < c,{n) < min{ i p'^in) ^ (i^g i)„-i^-„2-o(i) |_ 



If the initial diameter D is not 1 , then we must replace e hy e/D in the bounds for {n) . We easily 
check that the bound is essentially tight as long as e is not superexponentially small. Indeed, for any 
constant a > 0, there exist two constants b,c> such that, if e > p"", then (i) the communication 
count is at most (l/p)''"' and (ii) there exists an agreement system and a starting configuration for 
which the communication count exceeds [1/ pY^. Put more succinctly, 

Corollary 1.5. If e > p^'^''\ then Cein) = p^®("). 

Theorem 1.6. For any < e < p/n, an n-agent reversible agreement system e-converges to 
consensus in time 0{^n^ log 



Let 5 = maxj^t |A'j(t)| be the maximum degree of any node in the temporal network. The 
assignment qi = 5 \s valid and gives p = 1/5. The theorem implies e-convergence to consensus in 
0{5v? log ^) time. A similar result holds if the degree of any given node does not change over time 
and each agent moves to the mass center of its neighbors: Xi{t + 1) = {l/\Ni{t)\) YljeN^{t) Xj{t). If 
we now consider the case of a time-invariant graph, we retrieve the usual polynomial mixing time 
bound for lazy random walks in undirected graphs. 

The communication count is related to the total s-energy via the obvious inequality: 



In view of this relation, the two upper bounds in Theorem |1.4| follow directly from those in Theo 
simply set s 



rem 



1.1 



1 and s = n/ln^, respectively. Note that the second assignment can be 
assumed to satisfy s < 1, since it only concerns the case where 

ip-0{n) 

the bigger term in the right- 
1 / In ^ , and observe 



hand side of the expression in Theorem 1.4 For reversible systems, we set s 



that the number of steps witnessing a diameter in excess of e is at most e '^E^{s 



This bounds the time it takes for the diameter to dip below e and stay there forever (since it cannot 



grow); hence Theorem 1.6 



□ 



2 Applications 

We highlight the utility of the total s-energy by looking at five examples: opinion dynamics (^2.1 



social epistemology (^2.2); Kuramoto synchronization ( ^2.3[ ); bird flocking ( ^2.4[ ); and products of 
stochastic matrices ( ^2.5[ ). 
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2.1 Opinion Dynamics 



The Krause opinion dynamics model \W\ EI] is a sociological framework for tracking opinion polar- 
ization in a population. In its d-dimensional version, the bounded- confidence model, as it is often 
called, sets a parameter < r < 1 and, at time 0, specifies the opinions of n agents as n points in 
the unit cube [0, 1]*^. At time t > 0, each opinion x moves to the position given by averaging all 
the opinions that happen to fall within the Euclidean ball centered at x of radius r (or some other 
shape). Viewed as a multiagent agreement system, Gt consists of n nodes (the agents) with edges 
joining any two of them within distance r of each other. The dynamics is specified by 

* j€N,(t) 

where Ni(t) is the set of neighbors of node i in Gt, which as usual includes i itself. The system 



is known to converge [3 I27j . Theorem 1.4 allows us to bound how long it takes to reach 



equilibrium. Consider a Cartesian coordinate system. In view of ([5j[7]), we may set Pij{t) = l/|A'j(t)| 
and p = 1/n to make the opinion dynamics system along each coordinate axis conform to a one- 
dimensional multiagent agreement model ([3]). We can assume that the maximum diameter D along 
each axis is at most rn at time and hence thereafter. Indeed, by convexity, if along any coordinate 
axis the n opinions have diameter greater than rn, then they can be split into two subsets with no 
mutual interaction now and forever. Set £ = r/2 and let be the smallest t such that Gt consists 
only of edges in M'^ of length at most e. During the first dC^^^ (n) + 1 steps, it must be the case that. 



at some time t, the graph Gt contains only edges of length at most e. By Theorem 1.4 therefore 



Each connected component of Gt^ is a complete graph. To see why, observe that if opinion x is 
adjacent to y in Gt^ and the same is true of y and z, then x and z are at a distance at most 
2e = r, hence are connected and therefore at distance at most e at time ts- This "transitive closure" 
argument proves our claim. This implies that the opinions within any connected component end 
up at the same position at time + 1. Of course, when two opinions are joined together they can 
never get separated. The argument is now easy to complete. Either Gt^ consists entirely of isolated 
nodes, in which case the system is frozen in place, or it consists of complete subgraphs that collapse 
into single points. The number of distinct opinions decreases by at least one, so this process can be 
repeated at most n — 2 times. By this proves that Krause opinion dynamics converges in n'^('^^ 
time. We summarize our result. 



Theorem 2.1. Any initial configuration ofn opinions in the bounded- confidence Krause model with 
equal-weight averaging converges to a fixed configuration in n'^^") time. 

Martinez et al ^0] have established a polynomial bound for the one-dimensional case, d = 1. 
While extending their proof technique to higher dimension might be difficult, a polynomial bound 
could well hold for any constant d. We leave this as an interesting open problem. 



2.2 Truth-Seeking Systems 

In their pioneering work in computer-aided social epistemology, Hegselmann and Krause considered 
a variant of the bounded-confidence model that assumes a cognitive division of labor [17J. The idea 
is to take the previous model and fix one agent, the truth, while keeping the n — 1 others mobile. 
A truth seeker is a mobile agent joined to the truth in every Gt- All the other mobile agents are 
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ignorant, meaning that they never connect to the truth via an edge, although they might indirectly 
communicate with it via a path. Any two mobile agents are joined in Gt whenever their distance is 
less than r. (Using open balls simplifies the proofs a little.) Hegselmann and Krause [17] showed 
that, if all the mobile agents are truth seekers, they eventually reach consensus with the truth. 
Kurz and Rambau [22] proved that the presence of ignorant agents cannot prevent the truth seekers 
from converging toward the truth. The proof is quite technical and the authors leave open the 
higher-dimensional case. We generalize their results to any dimension and, as a bonus, bound the 
convergence rate. 

Theorem 2.2. Any initial configuration of n opinions in R"^ in the truth-seeking model converges, 
with all the truth seekers coalescing around the truth. If, in addition, we assume that the initial 
coordinates of each opinion as well as the radius r are encoded as 0{n)-bit rationals then, after 
^0{n) ^jj^g^ ^/jg ^jr-y^ffi seekers lie within a ball of radius 2""''" centered at the truth, for any 
arbitrarily large constant c > 0. Ignorant agents either lie in that ball or are frozen in place forever. 
This holds in any fixed dimension. 

Proof. Along each coordinate axis, a truth-seeking system falls within the fixed-agent agreement 



model and, as we saw in ^1.1, can be simulated by a (2re — l)-agent one-dimensional bidirectional 
agreement system with at most twice the initial diameter. (Note that the 2n — l agents do not form 
a truth-seeking system because there are no edges connecting the group of n original agents to its 



reflection.) Convergence follows from Fact 1.3, As we observed in the previous section, restricting 
ourselves to the equal- weight bounded confidence model allows us to set p = l/(2n — 1). (We could 
easily handle more general weights but this complicates the notation without adding anything of 
substance to the argument.) Kurz and Rambau [22] observed that the convergence rate cannot be 
bounded as a function of n and p alone because it also depends on the initial conditions (hence the 
need to bound the encoding length of the initial coordinates) . 

Set e = 2"''" for some large enough constant 6 > 0, and define tg as the smallest t such that 
Gt consists only of edges not longer than e. By the same projection argument we used in ([s]) and 
the observation that the initial diameter is 2^^"'\ t^ = nP^'^\ The subgraph of Gt^ induced by the 
mobile agents consists of disjoint complete subgraphs. Indeed, the transitive closure argument of 
the previous section shows that the distance between any two agents within the same connected 
component is at most 2e = 2^"*"" < r (the inequality following from the 0(n)-bit encoding of r), 
hence at most e. For similar reasons, the truth agent cannot join more than one of these complete 
subgraphs (referring here and below to the original system and not the duplicated version); therefore, 
all the subgraphs consist of ignorant agents, except for one of them, which contains all the truth 
seekers and to which the truth agent is joined. This truth group might contain some ignorant agents 
as well, ie, mobile agents not connected to the truth. For that reason, the truth group, in which we 
include the truth, is a connected subgraph that might not be complete. At time t^ + 1, the truth 
group has collapsed into either a single edge with the truth at one end or a coUinear 3-agent system 
consisting of the truth, a truth seeker, and an ignorant agent. (We refer to a single agent or truth 
seeker although it may be a collection of several of them collapsed into one.) All the other complete 



subgraphs collapse into all-ignorant single agents. By Theorem 2.1 there is a time 

to = te+ n^(") = (9) 

by which the all-ignorant agents will have converged into frozen positions unless they get to join 
with agents in the truth group at some point. 

Case I. Assume that the all- ignorant agents do not join with any agent in the truth group at any 
time t > t^: The truth group then behaves like a one-dimensional fixed-agent system with 2 or 3 
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agents embedded in 



We assume the latter, the other case being similar, only easier. We saw 



in il.l how such a system can be simulated by a one-dimensional 5-agent bidirectional system of at 
most twice the diameter. Recall that agents may represent the collapse of several of them, so we 
must keep the setting p = l/(2n — 1). The 5-agent system remains connected at all times (since 
its diameter cannot grow); therefore, by Fact 1.3 and Theorem 1.4, it /3-converges to consensus by 
(conservatively) time to + '^''^^^■'(log ^)^ time. By (9), this implies that, for any fixed cq > 0, the 
agents of the truth group are within distance 2~"'''° of the truth after n^^^^ time. 



Case II. Assume now that an all-ignorant agent z joins with an agent y of the truth group at time 
ti but not earlier in [te,ti). That means that the distance ||y(ti)2:(ti)||2 dips below r for the first 
time after t^. We want to show that ti < Iq + n'^^'^\ so we might as well assume that ti > Iq. Recall 
that tQ is an upper bound on the time by which the all-ignorant agents would converge if they never 
interacted again with the truth group past t^. Let L be the line along which the truth group evolves 
and let a be its (nonempty) intersection with the open ball of radius r centered at z{ti) = z{tQ). 
Note that a cannot be reduced to a single point. This implies that the shortest nonzero distance A 
between the truth and the two endpoints of a is well-defined. (By definition, if the truth sits at one 
endpoint, A is determined by the other one.) We claim that 

A>2-"°*"\ (10) 

Here is why. It is elementary to express A as a feasible value of a variable in a system of m linear and 
quadratic polynomials over m variables, where m is a constant (depending on d). The coefficients 
of the polynomials can be chosen to be integers over £ = n'^^"^ bits. (We postpone the explanation.) 
We need a standard root separation bound ^46j . Given a system of m integer-coefficient polynomials 
in m variables with a finite set of complex solution points, any nonzero coordinate has modulus at 
least 2"^')'°''"', where 7 — 1 is the maximum degree of any polynomial and i is the number of bits 
needed to represent any coefficient. This implies our claimed lower bound of 2""'°'"' on A. 

Why is £ = n*^^")? At any given time, consider the rationals describing the positions of the n 
agents and put them in a form with one common denominator. At time 0, each of the initial positions 
now requires 0{n?) bits (instead of just 0{n) bits). A single time step produces new rationals whose 
common denominator is at most n! times the previous one, while the numerators are sums of at 
most n previous numerators, each one multiplied by an integer at most n\. This means that, at time 
t, none of the numerators and denominators require more than 0(n^ + tn\ogn) bits. The system 
of equations expressing A can be formulated using integer coefficients with 0(n^ + tQulogn) bits, 
hence the bound of ^ = n'^^'^K Next, we distinguish between two cases. 



The truth is not an endpoint of a: Then there is a closed segment of L centered at the truth 



that lies either entirely outside of a or inside of it. By (10), the segment can be chosen of 



length at least 2~"''^*"'' . Setting cq large enough, as we saw earlier, the agents of the truth 
group are within distance 2""''°" of the truth after n*^^"^ time; therefore, ti < to + else 
the diameter of the truth group becomes too small to accommodate A. 

• The truth is an endpoint of a: Quite clearly, /3-convergence alone does not suffice to bound ti; 
so we reason as follows. When the truth group has /3-converged (for the previous value of /3), 
the only way its mobile agents avoided falling within a (in which case the previous bound on 
ti would hold) is if the truth group ended up separated from a by the truth (lest one of the 
mobile agents lay in a). By convexity, however, this property remains true from then on, and 
so z can never join y, which contradicts our assumption. 
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When agents y and z join in Gt at time t = ti, their common edge is of length at least r/3 
unless y on z has traveled a distance at least r/3 between and ti. In all cases, the system must 
expend 1-energy at least r/3 during that time interval. By Theorem this can happen at most 
n'^(")(3/r) = n*^*-") times. We can repeat the previous argument safely each time, even though the 
bit lengths will increase. At the completion of this process, we are back to Case I. □ 




Figure 6: Four coupled oscillators connected by four edges. 



2.3 Kuramoto Synchronization 

The Kuramoto model is a general framework for nonlinear coupled oscillators, with a dazzling array 
of applications: circadian neurons, chirping crickets, microwave oscillators, yeast cell suspensions, 
pacemaker cells in the heart, etc. Winfree's pioneering work on the subject led Kuramoto to formu- 
late the standard sync model for coupled oscillators \2>^\ I45j. The system consists of n oscillators: 
the i-th one has phase Oi and natural frequency oji. In its original formulation, the model is a 
mean-field approximation that assumes all-pair coupling. A more realistic assumption is to use a 
time-varying network to model communications. Considerable work has been done on this problem; 
see [131 EOl ESI |29l EH [351 ESI HH HZ] for a small sample. Further research introduced a time-1 dis- 
cretization of the continuous model |29l [33l 1371 SO] • Assuming all oscillators share the same natural 
frequency, a fixed phase shift yields the dynamics: 

where |-/Vi(t)| is the degree of i in the communication graph Gt-, which, as always, counts the self- 
loop at i (Figure [g]). As in [33], we also assume that all the agents' phases start in the same open 
half-circle. By shifting the origin, we express this condition as a — 7r/2 < 0i(O) < 7r/2, for some 
arbitrarily small positive constant a. This implies that 

sin(%(o) - em) = ai,(%(o) - em), 

for a/4 < Qij < 1. By ([5|; therefore, to make the dynamics conform to a bidirectional multiagent 
agreement system at time 0, it suffices to enforce the constraints: 

a 

Choosing p = ba/n for a small enough constant 5 > 0, we note that the constraints are roughly 
equivalent to < KAT < 1. By convexity, the angles at time 1 remain within [a — 7r/2,7r/2]; 
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therefore, our previous argument can be repeated to show that the synchronization dynamics fits 



within the bidirectional agreement model at all times. The result below follows from Corollary 1.5 
We note that it is impossible to bound the actual time to convergence unless we make assumptions 
about the temporal network. 

Theorem 2.3. Any Kuramoto synchronization system with n oscillators sharing the same natural 
frequency and initialized in an open half -circle e- converges after n*^*^") nontrivial steps, for any 
e > n"'^" and any constant c > 0. This holds regardless of the temporal communication network. 



2.4 Bird Flocking 

Beginning with Reynolds's pioneering work in the mid-eighties, bird flocking has generated an abun- 
dant literature, with a sudden flurry of interest in the last few years. Mathematically, flocking 
appears more complex than the previous agreement systems because the averaging and the commu- 
nications do not operate over precisely the same objects: it is the velocities that are averaged but 
the positions (ie, the integrals of the velocities) that determine the temporal network. Many models 
have been studied in the literature but most of them are variants of the following [6l [T2| [T9| HSj : 
given the initial conditions z{0) and ^(1), for any t > 0, 

( zit) = z{t - 1) + vit); 

\v{t+l)=P{t)v{t). ^ ' 

The vectors z{t),v{t) encode the positions and velocities of the n birds in M^: each coordinate of z{t) 
and v(t) is itself a three-dimensional vector. (These vectors are often expressed in M'^" via a tensor 
product; the notation here is easier as long as one remembers that the coordinates are themselves 
three-dimensional vectors.) The n-by-n stochastic matrix P{t) has nonzero diagonal entries and its 
other positive entries correspond to the edges of Gt', the communication graph Gt links any two 
birds within a fixed distance of each other. Intuitively, each bird averages out its own velocity with 
those of its neighbors in Gt'- all of its neighbors weigh equally in the average except perhaps for 
itself, ie, for fixed i, all nonzero pij{t)^s are equal, with the possible exception of Pii{t)\ all the entries 
in P{t) are rationals over O(logn) bits. 

It suffices to set p = n~^, for a large enough constant 6 > 0, to make flocking conform to 
the bidirectional multiagent agr eeme nt model, with v{t) encoding into a single vector the n points 



{xi{t), . . . ,Xn{t)). By Corollary 1.5, the system e-converges within n^^"^^ nontrivial steps for e > 
n~'^^ and any constant c > 0. We showed in [7] that the sequence Gt always converges to a fixed 
graph G, but that the number of steps to get there can be astronomical: it can be as high as a 
tower-of-twos of height on the order of logn, which, amazingly, is tight. 

Theorem 2.4. The velocities of n birds e-converge after n'^^") nontrivial steps, for any e > n"'^"' 
and any constant c > 0. The number of steps prior to the convergence of the temporal network to a 
fixed graph is no higher than a tower-of-twos of height O(logn); this bound is optimal in the worst 
case. 

2.5 Products of Stochastic Matrices 

Let V be the family of n-by-n stochastic matrices such that each P £ V satisfies the three standard 
constraints: (i) self-confidence (nonzero diagonal entries); (ii) mutual confidence (no pair Pij,Pji 
with exactly one 0); and (iii) nonvanishing confidence (positive entries at least p). Lorenz |27j and 



14 



Hendrickx and Blondel [18] independently proved the following counterintuitive result: in any finite 
product of matrices in V, each nonzero entry is at least /j'^^" ). What is surprising is that this lower 
bound is uniform, in that it is independent of the number of multiplicands in the product. We 
improve this lower bound to its optimal value: 



Theorem 2.5. LetV he the family ofn-by-n real stochastic matrices such that any P gV satisfies: 
each diagonal entry is nonzero; no pair pij , pji contains exactly one zero; and each positive entry is 
at least p. In any finite product of matrices in V, each nonzero entry is at least p^~^ . The bound is 
optimal. 



3 The Proofs 



It remains for us to prove: Theorem |1.1| (upper bound for s = 1 in ^3.1 upper bound for s < 1 
in ^3.2 and lower bounds in ^3.4), Theorem |2.5| ( ^3^ , Theorem 1.2 ( ^3.3[ ), and the lower bound of 
Theorem [L4l(pll). 



3.1 The General Case: s = l 



We prove the upper bound of Theorem 1.1 for s = 1. We show that En{l) < p ^^^^ by bounding 
the kinetic s-energy. 



Wingshift systems. We introduce a wingshift system, which provides a simpler framework for 
the proof. Since we focus on a single transition at a time, we write Oj, hi instead of Xi{t),Xi{t + 1) for 
notational convenience, and we relabel the agents so that < oi < • • • < a„ < 1. Given ai, . . . , a^, 
the agents move to their next positions hi,..., hn and then repeat this process endlessly in the 
manner described below. Let and r{i) be indices satisfying the following inequalities: 

RULE 1: 1 < i{i) < i < r{i) < n and {£ o r){i) < i < {r o i){i); 
RULE 2: a^(j) + 6i <hi < ar(.i) - 5i, where 6i = pia^f^i) - a^(j)). 

Each agent i picks an associate to its left (perhaps itself) and one to its right, £{i) and r{i), respec- 
tively. It then shifts anywhere in the interval [a^(j), a^^j)], though keeping away from the endpoints 
by a small distance 6i. This process is repeated forever, with each agent given a chance to change 
associates at every step. Any multiagent agreement system with parameter p can be modeled as a 
wingshift system: each agent picks its leftmost and rightmost neighbors as associates; note that the 
wingshift graph is sparser but now dependent on the embedding. Bidirectionality ensures rule 1: it 
says that the interval [£{i),r{i)] should contain i as well as all agents j pointing to ij^ By analogy, 
we define the total 1-energy of the wingshift system as ^ = Ylt>o where (with Oi denoting Xi{t)) 
Vt = ^^=1 {ar{i) ~ The desired upper bound En{l) = follows trivially from this bound 

on V: 



Theorem 3.1. The maximal total 1-energy of an n-agent wingshift system with unit initial diameter 
and parameter p is at most p~^^'^\ 



^This is necessary for convergence. Consider three agents ai = 0, 02 = |, as = 1, with ^(1) = r(l) = l{2) — 1 and 
r(2) = ^(3) = r(3) = 3. Agents 1 and 3 are stuck in place while agent 2 can move about freely forever. 
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Figure 7: A six-node wingshift system. 



As usual, we assume that p is smaller than a suitable constant. We need some notation to 
describe rightward paths in the wingshift system: r(i, 0) = i and r{i,k) = r{r{i,k — 1)) for k > 0. 
We define the distance between an agent and its right associate, = 0^(1) — a^. When traversing a 
rightward path i = r{i, 0), r(i, 1), . . . , r{i, k) etc, a sudden drop in Aj.(j j.) is of particular interest; so 
we introduce ri = r{i,ki), where 

ki = min | > | A^(j^fc) > j A^i^i^k+i) } ■ 

The agent r{ri) is called the stopper of i. Our interest in stoppers is that the nontrivial ones 
always move left, thus pointing to some obligatory motion in the system at the time step under 
consideration. 



Figure 8: The stopper of agent i is r{i, 3); rj = r{i, 2); and hi = 2. 



Lemma 3.2. The stopper u of any agent i satisfies au — 

Proof. By rule 1, ar(„) — a^(„) > A^. + A^. Since A„ < | A^., it follows by rule 2 that 

K < ar{u) - P{0.r{u) - < + A„ - p{/S.r^ + A„) 

< a„ + ((1 - p)p/2 - p)Ar, < a„ - (1 + p){p/2)'''+^ Ai . 

The last inequality follows from the fact that A,.. > (p/2)'^'Aj; it is not strict because both sides 
are equal if fc^ = 0. □ 

We bound V by tallying the kinetic 1-energy of the system, K = X^j>o Kt, where Kt = Y17=i 
Recall that and bi are the positions of agent i at times t and t + 1, respectively. 

Lemma 3.3. If K is finite, then V < n2"pi-"iC. 

Proof. Obviously, 

n 

t>0 i=l t,i t,i t,i 
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where = Oj — a£(i). By Lemma 



3.2 



E E^^ ^ E E (-) ' ("Kn) - ^(n)) < ^(-)"" E E 



\ai - bi\ . 



t>0 i=l t>0 i=l 

A mirror-image argument yields the same upper bound on Y2t i 



t>0 i=l 



□ 



The idea behind the proof. By symmetry, we can assume that at least half of the contribution 
to K is provided by rightward motions, ie, < ■{ bi — ai\bi > }. Thus we can conveniently 
ignore all leftward travel for accounting purposes. We use an "amortization" technique that involves 
assigning a credit account to each agent. Whenever the agent moves right it is required to pay for 
its travel cost by decreasing its account by an amount equal to the distance it travels. Credits are 
injected into the system only at time 0; if all travel is paid for and no account is overdrawn, then 
clearly the initial injection is an upper bound on ^K. The benefit of this approach is that accounts 
can borrow from one another, thus creating an "economy" of credits. The proof takes the form of 
an algorithm that drives the trading in a manner that keeps all accounts solvent; in other words, it 
is an algorithmic proof [8] . 



Agent i cannot, in a single step, move to the right by a distance greater than Aj. Lemma 3.2 



suggests a paying mechanism by which we charge its stopper u = r{ri) that travel cost; in other 
words, the leftward travel of u would pay for the rightward travel of the agents that claim u as a 
stopper. If u moves only to the left then its own travel distance is bounded by 1, and the charging 
scheme is essentially sound. But what if u zigzags left and right? The premise of charging u for the 
cost of i is that we know how to bound the cost of u. But, if u moves in both directions, we cannot 
bound its cost a priori (whereas we can if it only travels left). The solution is to look at u's own 
stopper u' and charge it. This may, in turn, force u' to charge u", etc. The "buck passing" evolves 
from left to right, so it must eventually stop. This picture suggests that agents should hold more 
credits the further to the right they are: indeed, our credit invariant will relate an agent's account 
to its rank. 



The algorithmic proof. At time t, the agents are ordered as < oi < • • • < a„ < 1. By using 
standard perturbation techniques, we can assume strict inequalities among all agent positions at 
all times. We maintain the following credit invariant: at the beginning of each time step, every 
agent i holds OjO* credits in its account, for some fixed parameter a, where again is shorthand 
for Xi{t). By way of illustration, consider the trivial case of two agents, one at and the other at 1, 
meeting at 1 — /? at the next step (we use ties for convenience). The system holds at time and 
(1 — p)(a + a^) at time 1. The difference exceeds the travel cost of 1 — /) if a is sufficiently larger 
than 1/ p. 

Our algorithmic proof involves setting up a simple data structure, a linked list, and moving 
credits around accordingly to specific rules. Let li be a stopper such that < a„. Consider the 
lowest-ranked agent h that claims u as its stopper. We build a doubly-linked list L„ consisting of 
u — h + 1 nodes, each one corresponding to an agent: h is at the head and u at the tail; scanning 
Lu takes us through the agents h,h + l, . . . ,u. The nodes scanned after v are called the antecedents 
of V. The rank s{v) is h plus the number of steps it takes to get from /i to f in L^. Ranks are implied 
by the list, so that inserting a node automatically adds one to the ranks of its antecedents. Liitially, 
the rank of the node v corresponding to agent k is just k and its position, denoted by a{v), is a^. 
The node following (resp. preceding) v in if it exists, is denoted by iiext(f ) (resp. prev(f )). We 
identify the node m with the highest-ranked agent such that a„i < b^- Let j3 = p/2 and a = 
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h 



m 



w 



step [1] 



u 



step [2] 



Figure 9: Credits are transferred via plain arrows in step [1] and dashed arrows in step [2]. Position 

a{v) = a/i, flm, &„, a„ for v — h,m, w, u. 



• For each stopper u from right to left, if 6„ < do: 

[1] Insert a new node w into right after m: next(w) -s— next (to) and 
next(TO) ^ w. Set a{w) = bu- For each antecedent v of w, transfer the 
account of v to prev(i;). Delete node u from L„. 

[2] For each node v from the new tail to next(/i): 

• transfer /3(au — a{v))a''^'"'^ credits from v to prev(t)); 

• keep {(3a{v) + (1 — /3)a„)Q;*^^^ credits in the account of v. 

[3] Move from a; to 6^ any agent i claiming u as its stopper, provided that 
bi > tti. Enforce all credit invariants. 

• Move from to bi any nonstopper i such that bi < a^. Enforce all credit 
invariants. 



Step [1] Since bu < a^, node u is the stopper of at least one node strictly to its left, so \Lu\ > 1 
and m is well defined. By assigning a{w) = bu, in effect we move the stopper u to its new position 
bu, right after agent m in L„. Shifting accounts one step backwards gets w to inherit the account of 
m + 1 and the new tail to receive the credits formerly at n. If m = u — 1, step [1] ends with the list 
in the same state as before except for a(tail); otherwise, the u — m — 1 antecedents of w see their 
ranks automatically incremented by one and, among them, the node for any agent k acquires the 
credit account of k + 1. (The alternative of keeping the list intact and shifting positions a{v) to the 
right works but breaks the immutable correspondence between nodes and agents.) To summarize, 
at the end of step [1], any node v in L„ ends up with a(v)a^^'"^ credits if v comes before w in the 
list and a{next{v))a^^'"^ otherwise|^ 

Step [2] We prove that all the credit allocations are feasible. Suppose that v is either w or an 
antecedent ofw. Agent v has a{'n.ext{v))a^^^^ credits. It receives /3(a„-a(next(t;)))a^('^^^*(^)) credits 
from next(i;) (which, by our notational convention, is zero if v is the new tail); it also gives away 

^By abuse of notation, a(next(v)) denotes if v is the new tail. 
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P{au — a{v))a^^^^ credits and keeps {/3a{v) + (1 — f3)au)(x^^^^ of them, for a total need of aua^^""^ ■ 
Since s(next(w)) = s{v) + 1, the transaction balances out if 

which holds because af3 > 1. Suppose now that v comes before w in the list. The only difference is 
that V now starts out with an account worth a{v)a^^^^ credits. The balance condition becomes 

a(z;)a^('') + /3(a„ - a(iiext(t.)))a^('')+i > aua'^^l 

which is equivalent to 

au — a{v) < a/3{au — a{iiext{v))). (12) 



To see why (12) holds, we turn to the wingshift condition — as, at some point, we must. Among the 
agents of Lu claiming u as a stopper and v as an antecedent, let z be the last one. (Note that z 
may not be equal to v, but because of h it is sure to exist.) Extending the notation in the obvious 
way, by Lemma 3.2 < (^)''^^Ha« - K)- If = 0, then 

au - a{v) . a« - flz 2 

^ -r~ ~ T' ^ ~ < ctP) 



On - a(next(t))) - 6„ au - bu p 



which proves (12). If A;^ > then, by the maximality of z, 

a(next(7;)) < 0^(2) < 0^(^(2)) < 



and condition (12) follows from 



au - a{v) ^ au - az _ + Ar(2) + - a^^riz)) ^ ^z + ^r(z) < ]^ _^ ^ 



a« - o(next(v)) Qu - ar(^z) ^^(2) + - ar(r(^)) \iz) 



Step [3] Having shown that all the accounts can afford the amounts specified in the second bullet 
of step [2], we now explain how they can pay for all rightward travel. We use an accounting trick, 
which is to move agents without crossing: all agents move continuously, one at a time. Should agent 
i bump into agent j, the latter completes the former's journey while i stops; the process is repeated 
at each collision. The main advantage of this scheme is that agents now keep their ranks at all times 
(of course they must also swap identities, thus becoming virtual agents). Unlike before, a rightward 
move Oj — )• bi will now entail the motion of all the virtual agents in the interval [oj , hi] and not only 
those claiming u as their stopper. By step [2], each node v is supplied with (/3a(w) + (1 — (3)au)ct^^'"^ 
credits, which is at least 

a{v)a'^'''> + (1 - p){au - a{v))a''^'"^ + a„ - a{v), 

because s{v) > 1. The three-part sum shows explicitly why virtual agent v, whose rank is now 
fixed, can move right by a distance of at least (1 — p){au — a{v)) while both maintaining its credit 
invariant and paying (comfortably) for the travel cost. Virtual agent v never needs to move further 
right than that. Why is that? The motion might be generated by v itself (if u is its stopper) or by 
the "push" from an agent on its left. Either way, at any instant during the continuous motion of 
virtual v, there is a causing agent i (perhaps v itself) whose corresponding interval [ai,6j] contains 
the position of v at that instant. Our claim follows then from rule 2. Indeed, 

bi < a,,(j) - p(ar(i) - a^(i)) < pa<?(i) + (1 - yo)a.r(i) < pa{v) + (1 - p)au < a{v) + (1 - p)(a„ - a{v)). 
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Returning all agents to their nonvirtual status, we observe that processing stopper u moves to 
the right only the agents that claim it as a stopper. Treating stoppers u in descending order from 
right to left means that none of the agents with u as their stopper has yet been moved (either to 
the left as stoppers or to the right) by the time we handle u. The last step in the boxed algorithm 
can only release credits — think of virtual agents to see why — and so, maintaining the corresponding 
invariants is immediate. This allows us to bound the kinetic 1-energy b}|^ 

1 " 

-K < ^Xi(0)Q^ < 2a" = 

j=i 



Theorem |3.1| follows now from Lemma 3.3 which completes the proof of the upper bound of Theo- 
rem [TTT] for s = 1. □ 



3.2 The General Case: s < 1 

We prove the upper bound of Theorem |1.1| for < s < 1 . We show that the total s-energy satisfies 
the recurrence: Ei{s) = and, for n > 2, 

En{s) < 2nEn-i{s) + (1 - {p/2ryEn{s) + n\ (13) 



We prove (13) "algorithmically" by describing a procedure to track the propagation of information 
across the temporal network and monitor its effect on the geometry of the system. All agents are 
initially dry, except for agent 1, which is wet. Every time a wet agent communicates with a dry 
one, the latter becomes wet. Once wet, an agent always remains so. Through the communication 
provided by the temporal network, water propagates from agent to agent. Bidirectionality ensures 
that, when the water ceases to spread to dry nodes, the interval spanned by the wet agents will 
be expected to have shrunken a little; in other words, communication acts as a spring that pulls 
recipients together. 



[1] Initially, all agents are dry except for agent 1. Set W{0) = {a;i(0)}. 
[2] For t = 0, 1, . . . ,cx): 

[2.1] Declare wet any agent adjacent to a wet agent in Gt- 

[2.2] W*{t) <~ W{t) U { positions at time t of dry agents just turned wet }. 

[2.3] Move each agent i from Xi{t) to Xi{t + 1). [ If no newly wet agent, 
then all motion within W{t) — W{t*) occurs in isolation from 
the n— |M^(i)| other agents. ] 

[2.4] W{t + 1) ^ { positions at time t + 1 of agents corresponding to W*{t) }. 



The set W{t) tracks the positions of the wet agents at time t. The auxiliary set W*{t) includes 
the positions at time t of the agents wet at time t + 1; it differs from W{t + 1) only in that the latter 

*A more sophisticated argument allows us to lower a to 0(l/p) and thus reduce the constant in the 0{n) exponent 
to 1; this sort of finotuning is not needed for the purposes of this paper. 
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gives the positions at time t + 1. Let ||P^(i)|| denote the length of the smallest interval enclosing 
W{t) and let {tk}k>i be the times t > 0, in chronological order, at which > |VF(t)| (ie, at 

least one dry agent turns wet at time Recall that p is smaller than a suitable constant. We 
show that: 



\\W{tk)\\<l-{^) . (14) 

The smallest interval [a,b] defining ||M^(tA;)|| is in [0, 1]. By symmetry, we can always assume that 
a + 6 > 1. Because = 0, we can also safely assume by induction that (14) holds up to tk] 



hence a > ]^{p/2)^. Since can increase only when at least one dry agent becomes wet, ie, at 

times of the form t = ti, we can prove (14) for tk+i by showing that + 1)|| < 1 — {p/2)^~^^. 

This easily follows from [0, ap) H W{tk + 1) = 0, so it suffices to prove the latter, which we do by 
contradiction. Consider an agent i contributing to + 1) with Xi^ti^. + 1) < ap. Agent i is wet 

at time tk + 1, so at least one agent in Ni{tk) was wet at time (possibly itself). This implies that 
Mi^t, > a and, by (§, 

Xiitk + 1) > (1 - p)mi,tk + pMi,tk > ap, 



which is impossible and proves (14) 



no 



ap 



a 



Figure 10: Bounding the interval spanned by wet agents. 

The set VF(tfc) can only gain agents, as k grows, but the set may stop growing before it absorbs 
all of them. When t is not of the form t^, the agents of TV(t) interact only among themselves, so the 
total s-energy expended during steps t^-i + 1, • • . , tfc — 1 is bounded by £'|iy(tj.)|(s) + (-s)- 
At time t = tk, the extra energy involved is 



^ |x.(t)-x,(t)r< (2 



Using obvious monotonicity properties, it follows that, up to the highest value of tk, the s-energy is 
bounded by 

n-l 

I / n \ I 



J2 \^Ei{s) + En_i{s) + (2) } ^ 



This includes the case where no tk exists. When it does and reaches its highest value t, if | VF(t+ 1) | < 
n then all the energy has been accounted for above. Otherwise, we must add the energy expended 



by the n agents past t. By (14), however, at time t + 1, the n agents fit within an interval of length 



l — {p/2)"'. By the scaling (power) law of the total s-energy, all we need to do is add {l — {p/2)^yEn{s) 



to the sum; hence (13). 



'Both W{t) and W*{t) are understood as multisets. Note that no might exist. 



21 



The case n = 2 is worthy of attention because it is easy to solve exactly. In the worst case, the 
two agents start at and 1 and move toward each other by the minimum allowed distance of p. 
This gives us the equation -^2(5) = (1 — 2/))'' £"2(5) + 1; hence, by (17), 

E2{s) 



1 1 

< 



1 - (1 - 2pY - 2sp 



(15) 



We now consider the case n > 2. By (17 13), 



En{s) < 



s{p/2Y 



By (15) and the monotonicity of En{s), we verify that the numerator is less than 3n^£'n-i(>s); 



therefore, for n > 2, by (15), 



EJs) < 



s{p/2y 



This proves the upper bound of Theorem |1.1| for s < 1 . 



□ 



Proof of Theorem 2.5, Recall that V is the family of n-by-n stochastic matrices such that any 
P £ V satisfies: each diagonal entry is nonzero; no pair Pij,Pji contains exactly one zero; and each 
positive entry is at least p. By ([s]), the entry of a product of t such matrices can be viewed 
as the position of agent i after t iterations of a bidirectional system with agreement parameter p, 
initialized with all the agents at 0, except for j positioned at Xj{0) = 1. Referring back to the boxed 
algorithm, we designate agent j as the one initially wet, with all the others dry. Let m(t) be the 
minimum value in W{t). At every time when W{t) grows in size, the minimum m(t) cannot 
approach closer than pm{t). Since \{tk}\ < n, either agent i stays dry forever and does not leave 
or it joins W{t) and cannot be smaller than minjm(t), which is at least p^~^. The lower bound 
proof suggests a trivial construction that achieves the very same bound and therefore proves its 



optimality. This completes the proof of Theorem 2.5 , □ 



3.3 The Reversible Case 



Our proof of Theorem 1.2 is based on a standard use of the Dirichlet form and classical spectral 
gap arguments [H [231 E2]- Let vTj = qi/^j Qj- We easily verify that vr = (tti, . . . , tt^) is the (time- 
invariant) stationary distribution of the stochastic matrix P = P{t) specified by (|4]): 



Pij 



i/qi 




if i = j; 
else. 



The argument focuses on a fixed step so we may drop t to simplify the notation. Let x = {xi, . . . , x„) 
and, for u,v ^ M", let {u,v)t^ = "^iTiuJvi. The dynamics is invariant under translation, so we may 
move the origin to ensure that {x, 1)^ = 0. Because vr is the stationary distribution, this property is 
time- invariant; in particular, {Px, l)7r = 0. Because P is reversible, we can decompose x = Yli^i'^^i 
in an eigenbasis {vi} for P orthonormal with respect to {■)-k] all the eigenvalues are real. Any positive 
Pij is at least l/qi > p. Let 1 = Ai > A2 > • • • > A„ > 2/9 — 1 be the eigenvalues of P, with the 
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labeling matching the fi's. Why the inequalities? Briefly, the gap is strict between the two largest 
eigenvalues because the graph is connected; the smallest eigenvalue is separated from —1 by at least 
2p because [P — pi) /{I — p) is itself a reversible Markov chain (with the same eigenvectors), hence 
with real- valued spectrum in [—1,1]. By Perron- Frobenius, if = max{A2, A^} then, by reversibility, 
T^iPij = T^jPji, and 

{X, X)^ - {PX, PX)^ = {X, (/ - P^)X)^ = '^^'^J{V^, U " P^)Vj)n = 0^(1 - A^) 

(16) 

> (1 - /x) = (1 - p)'Yaiaj{vi,Vj)n = (1 - p){x,x)t,. 

i i,j 

Because P is reversible and any nonzero nipij is at least p/n, it holds that, for any vector z, 



{z, {I - P)z)^ = T^iPijizi - Zjf > X] 



Set z = V2- By orthonormality, {z,z)t^ = 1 and (z, 1),^ = 0; therefore, z must contain a coordinate 
Za such that \za\ > 1 and another one, Zb, of opposite sign. Since Gt is connected, there is a simple 
path L connecting nodes a and b. By Cauchy-Schwarz, 

i-x, = {z,ii-p)z)^>£ E (^^-^^■)'^^( E \'^-' 

> {p/2n^){za- Zhf = p/2n\ 
Since A„ + 1 > 2p, it then follows that 



and, by ( |16[ ), 

{Px, Px)n < p{x, x)t, < (l - ^) 

Let E^{L,s) be the maximum value of the (diameter-based) total s-energy of an n-agent reversible 
agreement system such that (x, x)n = L at time 0. Since Gt is connected, qi > 2; therefore the 
diameter is at most 

2 max \xi\ < 2 IL/ minvrj < \/2Ln/ p ; 

i \l i 

therefore 

E^iL, s) < E^{{1 - pl2r?)L, s) + (2Ln/p)^/2. 

The total s-energy obeys the scaling law E^{aL,s) = a^/'^E^{L,s). The definition of E^{s) 
assumes unit initial diameter, which implies that {x^x)-,^ < 1, hence E^{s) < E^{l,s) and 

^ {2n/pYl^ 2n /2ny/2+i 

" ^""^ - 1 - (1 - p/2n2)-/2 - s\ p) 



which proves Theorem 1.2 This follows immediately from an inequality we shall use again later. 
For any < a, 6 < 1, 

(1 - a)^ < 1 - ah. (17) 



□ 
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3.4 The Lower Bounds 

We prove the lower bounds in Theorems 1.1 and 



1.4 



The case s < 1. We describe an algorithm An{a,b) that moves n agents initially within [a,b] 
toward a single point a + (6 — a)y{n) while producing a total s-energy equal to (6 — aYE{n,s). 
Clearly, £'(1, s) = 0, so assume n > 1. We specify A„(0, 1) as follows. Place n — 1 agents at position 
and one at position 1. The graph Go consists of a single edge between agent 1 at position 1 and 
agent 2 at position 0. At time 0, agent 2 moves to position p while agent 1 shifts to 1 — p. The 
n — 2 other agents stay put. Next, apply ^„_i(0,p) to the set of all agents except 1. By induction, 
we can assume that this brings them to position py{n — 1). Finally, apply An{py{n — 1), 1 — p) to 
all the agents. The operations of An leave the center of mass invariant, so if y{n) exists it must be 
1/n. Here is a formal argument. The attractor point y{n) satisfies the recurrence 

y{n) = py{n - 1) + (1 - py{n - 1) - p)y{n), 

where, for consistency, y{\) = 1. This implies that 

1 _ ^ 1 

y{n) y(n-l)' 

therefore y{n) = 1/n, as claimed. The total s-energy E{n, s) satisfies the relation: E{1, s) = 0; and, 
for n > 1, 

E{n, s) = p'E{n - 1, s) + (1 - py{n - 1) - pyE{n, s) + 1 
^ p'E{n-l,s) + 1 ^ pC'^-s)* 



1 - (1 - 2py - (1 - (1 - 2p)'*)"-i ■ 

Since p is small enough, {l — 2py > 1 — 3ps and E(n, s) > g^-"^ p-^i"-) ^ for any n large enough, s < sq, 
and fixed sq < 1- We observe that Algorithm An cannot start the second recursive call before the 
first one is finished, which literally takes forever. This technicality is easily handled, however, and 



we skip the discussion. This completes the proof of the lower bound of Theorem 1.1 for s < 1. □ 



The case s = 1. Suppose that each Gt consists of two nodes joined by an edge. The length of the 
edge can be made to shrink by a factor of 1 — 2p. We show that having n agents allows us to mimic 
the behavior of a 2-agent system with p replaced by (roughly) p": in other words, contraction can 
be made to slow down exponentially in n. Without loss of generality, we assume that n is an even 
integer 2m > 4. Our construction is symmetric by reflection along the X-axis about the origin, 
so we label the agents —m, . . . , —1, 1, . . . ,m from left to right, and restrict our discussion to the 
m agents with positive coordinates. (Equivalently, we could fix one agent.) The evolution of the 
system consists of phases denoted by ^ = 0,1,2, etc. At the beginning of phase 9, agent i lies at 
xi{0) = (1 - p'^f for i = 1 and af\ 

xi{6) = ^^-m+p'~H^-pn', 

for 2 < i < m. As usual, we assume that p > is small enough. The system includes a mirror 
image of this configuration about the origin at all times. Note that all the agents are comfortably 
confined to the interval [—2, 2], so the diameter D is at most 4. 

®We deviate slightly from our usual notation by letting the argument of Xi{9) refer to the phase of the construction 
and not the time t. 
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We now describe the motion at phase in chronological order, beginning with agent m. During 
phase 0, the first graph Gt {t = 6m) consists of exactly two edges: one joining m and m — 1 (with 
its mirror image across x = 0); the graph Gt+i joins m — 1 with m — 2 (and its mirror image); etc. 
The last graph in phase 6, Gt+m-i, follows a different pattern: it joins the two agents indexed 1 
and —1. Except for m, all of these agents (to right of the origin) are moved twice during phase 9: 
first to the right, then to the left. Specifically, agent 1 < i < m moves right at time t + m — i — 1 
and left at time t + m — i. We use barred symbols to denote the intermediate states, ie, the location 
after the rightward moves. At phase 6, 

^ . f x„^{e + 1) = a„x™„i(0) + (1 - a„)x„(0) = (1 - p"')xmi0); 

* ■ \ Xm-l{e) = ^Xm-m + lxm{e) = X„_i(0) + Ip"'-^! - p^' . 

We easily verify the identities above for am = {p — p"^~^^)/ {1 — p). For i = m — l,m — 2, . . . ,2, with 
Gt+ni-i joining agent i — 1 and i, the two moves are specified by: 



^ I x^{e + 1) = a,Xi_i(0) + (1 - ai)xi{e) = (1 - p-)xi(0) ; 

■ \ x,_i(0) = (1 - Pi)x,^i{e) + PiXi{9) = xi-i{e) + - p^f , 

where /3j = 1/(2 + p) and 

_ p 2(1 - p')/j™-*+i 
- 2T^ + (l-p)(2 + p) • 

Finally, at time t + m — 1, choosing ai = {p + 2p"^)/{A + 2p) allows us to write 

Gt+ra~i ■ xi{e + l) = -aixi{e) + (1 - ai)xi{e) = (1 - p'^)xi{9) . 

All the coefficients are 0(/>), so we can rescale p by a constant factor to make the dynamics 
conform to a standard one-dimensional bidirectional agreement system with parameter p (same 
with the diameter D). Obviously the system converges to consensus. In each phase 0, the union of 
the intervals formed by the edges of all of that phase's graphs Gt covers [—Xm{0)-, Xm{0)\] therefore, 
the total 1-energy is at least 



2Y^xm{e) = ^-^\^Y.^i-p-r>P 



■m 



This proves the lower bound of Theorem |1.1| for s = 1. For any positive e < 1/2, the length of the 
edge in Gt+m-i, which is 2xi{6), does not fall below e until 9 is on the order of p~'"log^, which 
establishes the lower bound of Theorem |1.4[ We note that the first agent oscillates around its initial 
position by roughly p/2 until 9 reaches so the kinetic 1-energy is, like the total 1-energy, 

exponential in n. □ 
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